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Abstract 

We investigate the conventional tight-binding model of L yr-electrons on a ring-shaped mol- 
ecule of L atoms with nearest neighbor hopping. The hopping amplitudes, t{w), depend on 
the atomic spacings, w, with an associated distortion energy V{w). A Hubbard type on-site 
interaction as well as nearest-neighbor repulsive potentials can also be included. We prove 
that when L = 4/c-t-2 the minimum energy E occurs either for equal spacing or for alternating 
spacings (dimerization); nothing more chaotic can occur. In particular this statement is true 
for the Peierls-Hubbard Hamiltonian which is the case of linear t{w) and quadratic V{w), 
i.e., t{w) = to — aw and V{w) = k{w — a)^, but our results hold for any choice of couplings 
or functions t{w) and V{w). When L = 4/c we prove that more chaotic minima can occur, as 
we show in an explicit example, but the alternating state is always asymptotically exact in 
the limit L ^ oo. Our analysis suggests three interesting conjectures about how dimerization 
stabilizes for large systems. We also treat the spin-Peierls problem and prove that nothing 
more chaotic than dimerization occurs for L = 4k + 2 and L = 4/c. 
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1. INTRODUCTION 

To derive the shape (and other properties) of molecules from first principles has been an 

actively pursued goal since the early days of quantum mechanics. Many of the insights into 
the structure of molecules like benzene and its relatives have, however, been obtained using 
drastically simplified models. The Schrodinger equation for all the nuclei and electrons in 
such a molecule involves many dozens or even hundreds of degrees of freedom and, therefore, 
simpler models with a reduced number of degrees of freedom are a necessity. In the case of 
benzene the introduction of the Hiickel model [1] played an important role. This model is 
standard textbook material for organic chemistry students (see e.g. [2]). London [3] used 
the Hiickel model to explain the large diamagnetic anisotropy of aromatic compounds and 
certain other materials quite successfully. A similar approach was used earlier by Jones in 
his work on bismuth and bismuth alloys [4] . 

In this paper we are interested in ring-shaped molecules of the type {CH)l, for even L, 
the so-called annulenes (sometimes called cyclic polyenes). The Hiickel model for [L]-annulene 
describes the L tt— electrons (one for each carbon atom) as hopping from one carbon atom 
to the next (tight- binding approximation). The carbon atoms are located at the L sites of 
a ring-shaped geometry, so L -|- 1 = 1. The Coulomb interaction between the electrons is 
ignored in the Hiickel model but we will include the Hubbard [5] on-site interaction as in the 
work of Pariser-Parr [6] and Pople [7] , in our study (a nearest neighbour repulsion can also 
be included). 

L L 

^^ = -E E h-c. +[/$^(n,T-i)(n,i-|) . (1.1) 

Here c]_,_]^^ and c^^ are the usual fermion creation and annihilation operators for a particle of 

spin a at site j. The number operators are nja — Cj^Cj^. The hopping matrix elements are 
real, but they can be multiplied by complex phases if a magnetic field is present. (In (1.1) we 
write (rija — |) instead of the customary njcr', the difference amounts only to the addition of 
trivial terms in the Hamiltonian, but these are convenient because they ensure — via hole- 
particle symmetry — that the expected particle number in the grand canonical ensemble is 
{N) — L, i.e., the "half-filled band".) The real parameter U represents the strength of the 
on-site interaction of two particles and the Hiickel approximation consists in setting U = 0. 
The set of tj's will be denoted by {tj}. The dependence of the Hamiltonian on the parameters 
will be made explicit, where needed, with a notation of the type H{{tj}), H{{tj},U) etc.. 
For a review of rigorous results on the Hubbard model see [8]. 

The hopping matrix element tj is a resonance integral between orbitals at the sites j 
and j + 1 and therefore depends on the real space distance Wj between these sites. A typical 
choice is tj = to — awj [9]. A more realistic choice would be tj = tQexp{—aWj) or some 
other rapidly decaying function of Wj. If one adds to the Hubbard Hamiltonian (even in the 
Hiickel approximation) a term describing the interaction between the carbon ions (the energy 
of lattice distortions) one obtains a model that, in spite of its extreme simplicity, serves very 
well to explain a number of phenomena in annulenes and linear {CH)^ (polyacetylene) [10]. 
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For a recent review see the references [11,12]. Although in more refined models the lattice 
distortions should be treated as quantum mechanical phonons [9], a more common choice, 
which we will follow here, is to describe them in the Born-Oppenheimer approximation by 
adding to the Hamiltonian (1.1) a classical potential of the form 

and to minimize, with respect to the Wj, the energy functional 

L 

where Xo{H) denotes the lowest eigenvalue of H. Often V is taken to be quadratic, but in 
this paper we will not assume that V is quadratic or that the dependence of the tj on the Wj 
is linear. 

The Peierls Instability [13] (discovered also by Frohlich [14] and seemingly independently 
by Longuet-Higgins and Salem [15]) states that for L sufficiently large and U not too large, 
the minimum will not be attained in a translation invariant configuration of Wj^s. By a 
straightforward computation (exact in the case U = or perturbative for small non-zero U) 
one can show that, for any fixed choice of the function V, there is an Lq such that for all 
L > Lq the ground state energy of H with 

Wj=wo + {-iy5 , (1.2) 

for small d, is lower than with the best choice of Wj — constant. In the context of {CH)x 
molecules this was discovered independently by Labhart [16] and Ooshika [17]. 

The phenomenon described by (1.2) is called dimerization in the physics literature. 
Unfortunately, this word has quite different connotations in the chemistry literature, but 
there does not seem to be a universally accepted terminology for (1.2) among chemists. The 
phrase bond- alternation [15] would be a more accurate description as far as chemists are 
concerned. We are obliged to make a choice here, and we shall use "dimerization" — in the 
hope that chemists will substitute "bond-alternation" for it in their minds. In fact we shall 
go further and declare a configuration to be dimerized, even if {1.2) holds with 5 = 0, i.e., the 
Wj's are translation invariant. In other words, a dimerized configuration is one with period- 
two translation invariance, and this includes period one as a special case. This convention, 
while a bit unusual, conveniently eliminates awkward locutions. 

For the Hiickel model {U = 0) and for L = 2 mod 4 (i.e., L = 6, 10, 14, . . .) the dimer- 
ization instability was shown very explicitly by Longuet-Higgins and Salem who also esti- 
mated the degree of dimerization 5 for realistic values of the parameters [15]. The physical 
mechanism for the occurence of lattice distortions is rather simple and quite universal for 
electron-lattice systems in one dimension. A lattice distortion of period l/^kp opens up a 
gap at the Fermi level, thus lowering the energy of the occupied levels. This was a basic 
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ingredient in Frohlich's theory of superconductivity. The calculations mentioned above show 
that for small distortions this lowering of the energy exceeds the positive contribution to the 
total energy of the lattice distortion itself. In the infinite volume limit the elastic energy per 
bond is quadratic in the parameter 6 of (1.2), while the electronic energy per bond decreases 
by an amount proportional to 5^ log 5 for small 5. As far as we know the first author who 
mentions this logarithmic behavior is Prohlich in [14]. 

A common feature of the works mentioned above is that they always show that certain 
instabilities exist. The question of the "stability of the instabilities" is usually not raised, 
much less resolved. The main purpose of this work is to study the stability and instability of 
lattice distortions and to obtain rigorous statements about the true energy minimizing state. 
Longuet-Higgins and Salem [15] raised the question of the possible occurrence of higher 
periodicities but expected them to be unimportant. We prove that in the case L = 2 mod 4 
(an assumption they made for technical reasons) they were right: indeed nothing else than 
periodicity two occurs for L = 2 mod 4, but when L — mod 4 other instabilities may, in 
fact, occur. 

For the Hiickel model with t = to — aw and V{w) — k{w — a)'^, Kennedy and Lieb 
[18] showed that the true minimizing configuration is always of periodicity two, i.e., wj is of 
the form (1.2). Moreover, the minimizing configuration is unique (up to translations). It is 
remarkable that this result holds for all even L, in contrast with what will be proved here 
for general functions V. 

We note parenthetically that the models we study are also used to describe electrons in 
mesoscopic metallic rings [19] . These rings are typically two orders of magnitude larger than 
the largest annulenes produced in the laboratory so far, and therefore the magnetic fields 
needed to obtain a flux through the ring of the order of a flux quantum are much smaller 
and experimentally accessible. The main issue is to calculate the persistent currents in such 
a ring threaded by a magnetic flux. Depending on the electron number, among other things, 
these currents can be paramagnetic or diamagnetic. Although most authors assume uniform 
hoppings t for this problem, there is an intimate connection between these persistent currents 
and the stability of dimerized configurations of Vs. 

In order to proceed to a precise formulation of our results it is first necessary to confront 
some possibly confusing questions about the phases of the t^'s. Let us write 

tj = \tj\exp[iej] 

with — TT < dj < TT. We note, first, that by a simple unitary gauge transformation of the type 
Cj^ — > e'xp[i(pj]cj^ and cj^ exp[— zc/jjjcj^, the operators nja- and the energy levels of H are 
unchanged, but the tj's change to tj tj exp[i{(pj — (pjj^i)]. Consequently, the energy levels 
depend only on the "total fiux", $, of the t/s defined by 

L 

$ = argument (JJ^ tj) , 

with — TT < $ < TT. 
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The next consideration comes from the physics of our model. The t^'s are real, unless 
there is a superimposed magnetic field, in the absence of which each 9j = or tt, depending 
on whether tj is positive or negative. Both signs can occur. Furthermore, it is understood 
that the t/s remain close to some nonzero value in the physically realistic situation. Without 
loss of generality we can choose this value to be positive. With this in mind, we shall assume 
that the elastic energy, as a function of t, defined by 

f{\t\)^mi{V{w)\t{w)=t} , 

depends only on \t\ and not on the sign of t. (We will impose below in (1.6) some other very 
mild physical conditions on /.) This is not to say that negative Vs will not be permitted. 
They will — as wc explain next. 

Our procedure will be to study the total energy 

£Lm)^Xo{H{{tj})) + J2fi\tj\) , (1-3) 

j=l 

and to evaluate its minimum 

EL = mm£L{{tj}) , (1.4) 

where the minimum is taken over all choices of {\tj\} as well as the phases {Oj}. Thus, 
we allow negative tj's and give them the same elastic energy as \tj\. Since the tj's are not 
expected to pass through zero in the physical case, this is no real restriction. In terms of the 
tj, dimerization means that the configurations minimzing SL{{tj}) are of the form 

tj = to + i-iy5 . (1.5) 



It will turn out that for all L our theorems show that the energy minimzing tj's are real 
(modulo a gauge transformation). But for L = mod 4 an odd number of them must be 
negative for an energy minimizer — thereby precluding a dimerized state. If, on the other 
hand, we wish to attribute the occurrence of an odd number of negative t^ 's to the presence 
of an external magnetic field then we can say that a dimerized state minimizes the energy 
in the presence of a field of fiux tt. Given that no such field is really present, we will have 
to conclude that the energy mimimum will not be a dimerized state when L = mod 4. In 
this case we can redefine our problem by restricting the minimum in (1.4) to positive tj^s. 
Do we get a dimerized state then? We do not know the general answer but we have some 
conjectures about this question. If f(t) = k{t — a)^, as in [18], we know that a dimerized 
state is, indeed, the minimum. 

The physical conditions we impose on / are 

i) f is continuous. 

' (1.6) 

ii) f{t) > Ct, for large t and for some constant C > 4. 
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The continuity of / is not really necessary for the validity of our results. In fact in the 
Appendix (Lemma 13) we show that / can always be replaced with a convex (and hence con- 
tinuous) function without changing the minimizing configurations. Property ii) is necessary 
to have a stable minimum at all, e.g., to prevent the molecule from collapsing to a point. 

We minimize the ground state energy £l not only with respect to the parameters tj but 
also with respect to the number of electrons. In Section 2 we prove that if [/ 7^ and all 
tj are real, the ground states of (1.1) will in fact necessarily have N = L. In this case all 
results below are equally valid for minimization of the energy with the particle number fixed 
at half-filling. The condition [/ ^ is necessary in the case of rings of size L = mod 4, for 
with t/ = the model with translation invariant t's has ground states with particle number 
ranging from L = — 2 to L -|- 2, due to zero eigenvalues of the one-particle Hamiltonian. 

Relation with the Flux-Phase Problem. 

The minimization with respect to the flux $ of ££,({tj}), for fixed \tj\, is a generalization 
of the so-called fiux-phase problem [20,21] for rings. In the case C/ = (as in the original 
formulation of this problem), but arbitrary, fixed, positive {tj}, this problem was solved by 
Lieb and Loss [22]. The arguments of this paper provide a way of solving the problem for 
rings and all values of U (Corollary 8) . For rings of length L = 2 mod 4 the optimal (= energy 
minimizing) flux is 0. When L = mod 4, the optimal flux is <^ = tt. This can be extended 
to higher dimensions [23] when some geometric periodicity is present. In particular, [23] 
proves the conjecture that the optimal fliix for the Hubbard model on the two-dimenisonal 
square lattice is tt through each plaquette. 

Main Results. 

The following theorem, which we prove in Section 3, says that, when L = 2 mod 4 and the 
flux is zero and when L — mod 4 and the flirx is tt, the Peierls instability is itself stable in 
the sense that no other lattice distortions than period two can occur. The theorem does not 
state that distortion always happens. Depending on the parameters in the Hamiltonian, the 
function / in particular, the energy minimzing configuration may or may not be translation 
invariant. In real molecules both situations occur. Benzene, e.g., has all C — C bond lengths 
equal (within the precision of today's measurements and in agreement with the analysis of 
Labhart [16] and Ooshika [24]), while for large rings or long chains the Peierls instability 
will necessarily lead to two different bond lengths. 

1 Theorem (Minimizing configurations). 

i) If L = 2 mod 4, the minimum of £l is attained in a dimerized configuration of the form 
(1.5) in which the tj 's are all positive (or all negative). 

ii) If L = mod 4, the minimum of £l is attained configuration in which the \tj\ 's satisfy 
(1.5) but the flux $ equals n. 

That the cases L — 2 mod 4 and L — mod 4 behave differently with respect to the fiux 
was noted a long time ago (see e.g. the "Hiickel rule" in [2].) The U = one-electron problem 
has zero-energy states when L = mod 4 but not when L = 2 mod 4. In the first case this 



6 



leads to diamagnetic, in the second case to paramagnetic response to magnetic fields, as was 
noted in [3] and as is observed experimentally in the annulenes [25] . 

2 Theorem (Uniqueness) . // there is only one ( up to translation and gauge transforma- 
tions) dimerized configuration satisfying the conclusion of Theorem 1, then there are no other 
energy minimizing configurations. 

We define three additional minimization problems for the L = mod 4 case. 

• Zero- flux problem: We restrict the minimization in (1.4) to those {tj} with zero fiux, 
i.e., $ = 0. 

• Period-two problem: In addition to $ = we make the further restriction that the tj's 
must be dimerized. 

• The infinite- chain problem: Let y) denote the dimerized configuration with ^2^+1 = x 
and t2j = y, and minimize the energy per site defined by 

e{x,y)= lim ySLid{x,y)) . (1.7) 

L— >oo L/ 

The minimum energy in the first case is denoted by El, and the t^'s by The second 
energy is denoted by E^^^^^^^ and the t/s by t°dimer- Clearly 

-'^L.dimer ^ 

and one of our conjectures below concerns the case of equality. The following, however, 
shows that the three problems become asymptotically the same (to within two powers of 
L~^, instead of merely one) as L — > 00. 

3 Theorem (Asymptotic dimerization). 

i) Assume that there is a unique energy minimizing dimerized configuration d{x,y), x > y, 
for the infinite- chain problem. Then any sequence of minimizing configurations {t^} of 
the zero-flux problem converges to a dimerized configuration as L — > 00. Namely, 
lim max|t^ ~^?dimerl ~ max|t^ — d{x,y)j\ = , 

where we have adopted the convention that ti is the largest of the tj 's for all configura- 
tions. 

a) For any even L, the energies satisfy 

0<Et 

4imer ~ ^£ — constant X L ^ 

In section 5 we discuss the situation of finite rings of length L = mod 4. Instabilities 
other than period-two can exist for arbitrarily large rings if the function / is chosen appro- 
priately. An optimally dimerized ring of size L — mod 4 can have other instabilities and it 
seems that the strongest of those actually has wavelength = L, and is not, as one could have 
expected, of period 4 or some other higher periodicity. We illustrate this with an explicit 
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example of a ring of 8 sites for which the minimum is not period 2. A perturbative analysis 

indicates that the conditions on / for a further instability to occur become more stringent 
as the system becomes larger. In particular the dimerization itself has to be smaller as L 
becomes bigger. This leads us to the following two conjectures: 

Conjecture 1. For every L — mod 4 there exists a function / (depending on L), satisfying 
conditions i) and ii) in (1.6), such that the energy minimizing configurations for the zero-flux 
problem are not period two. 

Conjecture 2. For any fixed / satisfying (1.6), and such that there is a unique minimum 
for the infinite-chain problem, there is a critical size Lc such that for all even L > Lc all 
zero-fiux minima are period two. 

The Spin-Peierls Problem. 

Another closely related question is the spin-Peierls problem. For U large, and half-filling, 
the Hubbard model reduces to the Heisenberg antiferromagnetic spin chain in second order 
perturbation theory [26]. The spin chain also has an instability when coupled to lattice 
distortions, as was pointed out by Chesnut [27] and by Beni and Pincus [28,29]. More 
explicitly one looks for the configurations of coupling constants that minimize the smallest 
eigenvalue of the Hamiltonian: 

<i,j> <i,j> 

where the sum is over nearest neighbour pairs in a finite subset A of the hypercubic lattice, 
and f{Jij) is the elastic energy as before. In one dimension Cross and Fisher [30] computed 
the exponent governing the ground state energy density e{S) of the system with alternating 
couplings: 

j^ = i + {-iy5 , 

with Jj = Jjj+i- They found that, up to logarithmic corrections, 

|e(5)-e(0)| _ 

Numerical confirmation of this exponent and evidence for a logarithmic correction was re- 
ported in [31]. The important point here is that 4/3 < 2, which implies the instability of 
the translation invariant configurations under period 2 perturbations, provided / has a finite 
second derivative. It was recently proved that this exponent is identical to a critical exponent 
of the two dimensional 4-state Potts model [32]. 

The question addressed in this paper is again to determine the nature of the minimizing 
configurations. It turns out that the true minimum is dimerized for any (even) system size and 
in any dimension. The 2 mod 4 versus mod 4 dichotomy does not arise for the Heisenberg 
model! 

4 Theorem (Dimerization for spin-Peierls). Let A C Z'^ 6e a rectangular box of even 
size in all coordinate directions and with periodic boundary conditions. Then there is an 
energy minimizing configuration of Jij 's which is of periodicity 2 in all coordinate directions. 
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This theorem prompts one more conjecture. 

Conjecture 3. Suppose that a minimizer for the zero-flux problem for some / and L = 
mod 4 and some C/q ^ is a dimerized state. Then for every U > Uq there is a dimerized 
minimizer for the zero-flux problem with the same / and L. In particular, the truth of this 
conjecture, when combined with the U — result in [18], would imply that, with a quadratic 
/, dimerization always occurs, even in the L = mod 4 case. 

Extensions of the results presented in this paper. 

1) It has been argued that the Hubbard on-site repulsion does not describe the Coulomb 
interaction between the electrons accurately enough (for a discussion see e.g. the reviews 
on this topic in [33] and references therein). While the Hubbard term alone, when not too 
large, seems to have the effect of enhancing the dimerization [34] , taking into account some 
nearest neighbour interaction terms Kivelson, Su, Schriefer, and Heeger [35] found that the 
Coulomb repulsion in fact suppresses the dimerization in polyacetylene. Without entering 
into this discussion here we would like to point out that without substantial modification 
one can prove the results of Theorem 1 and Theorem 2 for a Hamiltonian that includes a 
repulsive nearest-neighbour interaction of the form 

3 

The — I's are inserted to keep the chemical potential tuned at half-filling. We require W{t) > 
0, but no assumption on the dependence of on t is needed. Theorem 1 holds unchanged, 
but we can extend the proof that the grand canonical ground states are necessarily half-fllled 
only under the additional condition that W{tj-i) + W{tj) < Uj, for all t's in the relevant 
range. Nearest neighbour exchange terms and certain longer range interactions can also be 
treated. 

2) Instead of a homogeneous on-site repulsion it is sometimes more realistic to have 
different U for different (classes of) sites, e.g. in a periodic manner. Theorem 1 extends to 
such situations in the following sense: the minimum energy will be attained in a conflguration 
of tj which is invariant under all reflections through planes that intersect two oppossite bonds 
on the ring (as indicated by a dashed line in Figure 1) and that leave the configuration of the 
Ui invariant. As an example consider [18]-annulene. Of the 18 hydrogens in this molecule 
6 reside inside the ring and 12 outside, according to the pattern: 1 inside, 2 outside, 1 
inside, etc. A possible way of taking the effect of this periodic configuration of H-atoms 
into account, is to add a periodic one-particle potential to the Hamiltonian. In the present 
example this potential would take two different values and depend on the site according to the 
pattern: (vi) — (f i, f i, f2, 'f^i, i^i, • • •)• general result then implies an energy minimizing 
configuration of C - C bonds of the form ABC AC B ABC AC B The observed bond 
lengths in [18]-annulene indeed satisfy this pattern (they are A = 1.419A, B = 1.382A,C = B 
[25].) 

3) Hubbard models with spin-dependent hoppings, i.e., tj = \tj\mj,mj G SU(2), have 
been considered in the literature in order to study the effects of spin-orbit coupling [36,37,38]. 
By virtue of the SU(2) gauge symmetry of the Hamiltonian one can diagonalize the spin de- 
pendence of the hoppings up to an SU(2) flux, which is deflned up to an SU(2) transformation 
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by the matrix M = mim2m^ ■ ■ -rriL. The methods of this paper can be used to prove the 
existence of a minimizer having dimerization and M = I when L = 2 mod 4 and when 
L = mod 4 a minimizer has M = — I and dimerization. This is completely analogous to 
the U(l) case treated in Theorem 1. 

4) Mattis and Langer [39] introduced a simple model to study the interaction between 
the electrons and a phonon instability of the Peierls type as a function of the temperature. 
Our methods can also be used to prove that in this model the electron state has periodicity 
two and that its correlations indicate dimerization. For a rigorous study of the Kohn anomaly, 
which accompanies the phase transition in the Mattis-Langer model, see [40]. 

Generalizations that are not considered in this work include the interaction between 
the electrons and lattice in excited states, the effects of doping, and higher dimensional 
models. Low-lying electronic excited states can induce further (e.g. soliton-like) distortions 
of the lattice [9,41,42]. In doped polyacctylcne an interesting semiconductor-metal transition 
occurs [43] . Away from half-filling, instabilities other than period-two will naturally develop 
(see e.g. [44] and references therein). In two dimensions, e.g., on a square lattice, 2kF 
instabilities may develop either in the coordinate directions or in the diagonal directions 
[45,46]. This may be related to the non-period-two instabilities that can occur in rings 
whose size is a multiple of 4, e.g., in the elementary plaqucttes of the square lattice. Another 
new feature of higher dimensions is that the breaking of translation invariance can occur at 
non-zero temperature. This has been rigorously shown to occur in the Falicov-Kimball model 
[47,48] and in the Holstein model [49] . 



2. THE GRAND CANONICAL HUBBARD MODEL AT HALF-FILLING 



In the definition of the energy functional (1.3) we used the lowest eigenvalue of the grand 
canonical Hubbard Hamiltonian at half-filling. Our goal in this section is to prove that the 
corresponding ground state necessarily has = L if we assume the tj's are real (as will later 
be shown to be the case). It is then obvious that minimization of the energy with respect 
to the parameters tj and the number of electrons, produces exactly the same minimizing 
configurations as the restricted mimimization with the particle number fixed at A = L. The 
arguments in this section are not restricted to one dimension. We believe they are interesting 
in their own right and can be of interest in a more general context. Therefore, we temporarily 
consider a more general setup on a general finite lattice A, the number of whose lattice sites 
is denoted by |A|. 

Consider the grand canonical Hubbard Hamiltonian (1.1) for spin 1/2 fermions on A, 
and with arbitrary chemical potential /i: 



The tij are asssumed (in this section only) to be real and symmetric, i.e tij = tji. A is 
assumed to be connected, i.e., given i and j in A we can find a sequence i — ii,i2, ■ ■ ■ ,im — j 
such that tij^, 12^2, is ■ ■ ■'ti^-i,im 7^ 0- The Uj and /j, are also real parameters. As before the 
dependence of H on its parameters will be made explicit with a notation of the type iif^. 



i,jeAa=-\,l 
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H{{tij}), Hi^{{tij}, {Uj}) etc. When H appears without subscript the chemical potential is 
= 1/2, corresponding to half-filling. 

The spin-up and spin-down particle number operators are defined by 

j&A 

where a =|, | and n^v = cj^c^^. Note that and Ni each commute with if, and hence 
their eigenvalues are good quantum numbers. H also commutes with the total-spin operators 
(which generate the global spin rotations) given by 

Therefore, by applying S'^ or 5"" an appropriate number of times, one can transform any 
ground state of H into a ground state with 

iV,-7V,=(° if iV is even . . 

^ ^ l±l ifiVisodd ^ ^ 

We say that A and t are bipartite if the sites of A can be written as the disjoint union of 
two subsets; i.e., A = A^ U Ab, in such a way that the matrix t of hopping matrix elements 
respects this structure, namely tij = if i G A^ and j G A^ or if i G A^ and j E As- In this 
case the Hamiltonians Hi/2{{tij}, {Uj}) and Hi/2{{tij}, {—Uj}) are unitarily equivalent, i.e 
at half- filling (/i = 1/2) a global sign change of the potential does not affect the spectrum. 
A unitary transformation implementing this equivalence is the particle-hole transformation 
mapping 

S'i ' Si' Si ' Si 

where e{j) = 1 for j G A a and e{j) = —1 for j e Ab- By the same transformation the spin 
operators are mapped into the so-called pseudospin operators defined by 

s"- = E^(.^')s-TSi' ^' = E^(-^>JisV = -N^- 

These operators therefore also commute with the Hamiltonian, and they also commute with 
the spin operators and generate an additional SU(2) symmetry. Under the particle-hole 
transformation N-^ and are transformed into |A| — N-^ and respectively. It is then 
obvious that by applying S"^ to any ground state we can obtain a ground state with 

^T + A^i = l!tLT (2-2) 
' I A ± 1 if A/^ - A IS odd ^ ^ 
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As the spin and pseudospin operators commute (2.1) and (2.2) can be reahzed simulta- 
neously. For an even bipartite lattice and ^ = 1/2 this implies that there is a ground state 
of H with Ni; = Ni = |A|/2 or A^^ = ± 1 = |A|/2. 

From the next lemma it will follow that if in addition all Uj are nonvanishing and of the 
same sign, all ground states have even particle number and satisfy N-^ = — |A|/2. Note 
that the Lemma itself does not require A to be bipartite. 

5 Lemma. // all tij are real and Uj < for allj, then the ground state space of H{{tij}, {Uj}) 
contains a state with — N^. If Uj < for all j , all groundstates satisfy N-^ — ^l- 
particular the total number of particles in the ground state is even and the total spin is zero. 

Proof : We can introduce new operators ^ defined by 

H has then the same form in terms of the c's as in terms of the c's, but now the cf , commute 
with the cf-f. 

Our Hamiltonian is now an operator on H,^ ® Hi, and can be written in the form 

H = T0l + l0T + J2Uj{nj-i^)(S){nj- jj) 

where T = — Y2i jeh ^ij^l^i- (Because spin f and spin J, operators act identically on different 
tensor factors in the Hilbert space, we can omit the spin index.) T and the Uj are real in the 
canonical basis of localized particles. Therefore we can apply Lemma 14 of the appendix and 
the discussion thereafter to conclude that H has a ground state with A| = N^. 

If Uj < for all z, there cannot be a ground state with ^ N^. Indeed the last statement 
of Lemma 14 implies that for any ground state O 

Uj{^ I {nj,i - /x)n) = Uj{n I {nj,i - n)n) 
If Uj < this implies = N^. Q.E.D. 



6 Lemma. // all tij are real and bipartite, |A^| = \Ab\, and either all Uj > or all Uj < 0. 
Then the ground state of H with fx = 1/2 is unique and has = = |A|/2. 

Proof : It is sufficient to consider the case of all Uj > 0. The case Uj < then follows because 
the particle- hole transformation changes the sign of the Uj, and the properties = |A|/2 
and Ni = |A|/2 are unchanged under the transformation. Let us therefore consider the case 
Uj > 0, for aU j. 

We first apply Lemma 5 to the particle-hole transformed Hamiltonian which has all Uj < 0. 
This tells us that = |A| — in all ground states. In particular the total number of par- 
ticles is |A| which, by assumption, is even. Then, by Theorem 2 of [50] (the grand-canonical 
version at half-filling does not need the homogeneity of the potential) and the assumptions 
on the lattice, N-^ = = |A|/2. [Note: Theorem 2 of [50] requires Uj = constant. That is 
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the case when the interaction is Tij-iUji as in [50]. If we use the formulation of (1.1), however, 
with (n^i- — l/2){nji — 1/2), the same proof works without the need of constant Uj.] Q.E.D. 

For translation invariant tij on regular lattices the case Uj = is of course exactly 
solvable. It is rather straightforward in such cases to determine the ground state degeneracy 
"by hand". The ground state is unique if and only if the eigenvalues of the matrix (Uj) are 
all non-zero. 

It is straightforward to extend Lemma 5 to Hamiltonians that include interactions of the 

form 

a 

for arbitrary real functions F and G of the local particle numbers. One can use this to prove 
a slightly weaker form of Lemma 6 for a class of models with an additional term of the form 

added to the Hamiltonian H. For these models our methods also prove that 

i) if Wij > and {Uj — ^^Wij) > 0, for all j, then there is a ground state with 
N^=N^ = |A|/2. 

n) if Wij > and {Uj - Wij) > 0, aU ground states have = = |A|/2. 



3. GENERAL RESULTS BASED ON REFLECTION POSITIVITY; 
PROOFS OF THEOREMS 1 AND 4 

We now return to the half-filled Hubbard model (A^ = |A|) on a ring of even length |A| = L. 
It will be convenient to express the dependence of the energy functional on {tj} in a more 
explicit way. We write 

£Ut^^\t^"^\&^) = Xo{H{{tj})) + f2fm > 

J=l 

where t^^\t^^\ and t^'^^ denote a partition of the complex tj's (with t* denoting the complex 
conjugate of tj) into three groups as follows (also see Figure 1): 

= {h,...,tL/2-l), t^"^^ = {tL,tL/2), t^'^ = {tL-l,...,tL/2+l) ■ 

The reflection of a configuration {t^'-\t^'^\t^'^^) through the plane intersecting the bonds 
{L, 1} and {L/2,L/2 + 1} is the configuration {&\t^'^\t^^^), i.e., with t^^) ^ &\ 

Throughout this section we assume a constant potential Uj = U. The following lemma 
can be formulated to include non-constant potentials but this is not needed for our purposes. 
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7 Lemma. Assume a constant potential Uj = U. If L = 2 mod 4, assume that t]^/2 and tL 
are real and nonnegative. If L = mod 4, assume that ^ and < 0. The other tj 
are allowed to be arbitrary complex numbers. Then the energy functional satisfies 

SL{t^^\t^'^\t^'-^) > ^ (SL{t^^\t^'^\t^^>)+SL{t^''>,t^'^\t^''^)j . (3.1) 

Proof : As in the proof of Lemma 5 we consider the Hamiltonian (1.1) as acting on Ti^ ^T^i, 
which amounts to considering operators referring to opposite spins as commuting. For a =t, i 
define a set of "spin" operators using a Jordan- Wigner transformation as foUows: 
c3 _ t „ _ 1 c . _ 9 c3 

With these definition and the anticommutation relations of the ^ and ct ^ it is straigth- 
forward to check that ^, S^^, and 5""^ satisfy the standard SU(2) commutation relations 
and that they commute for difi'erent indices i,a. In terms of these operators, and using the 
commutation relations, the Hamitonian can be written as 

tL E '5ita(Si,c.---SL,.)^Z,.+ h.c. . 

By performing a rotation by tt about the 2-axis at every other site and for both a and 
o^ =1 we find that this spin Hamiltonian is unitarily equivalent to: 

L-l L 

^ = E E h.c. +E^^-4T^a 

+ {-if^'tL E '5t.(Sl,. • • • SL/2,a)StA^L,a ■ ■ ■ SL/2+l,a) + h.C. . 

Under the conditions stated in the proposition tL/2 > and {—l)^^^tL < 0. We can now 
apply Lemma 14 of the appendix to H with 

L/2-1 L/2 

^ = - E E h.c. +uj2shsli , 

^ = - E E t3S|+^,.Sl + h.c. +U f2 ^hS^l ' 

i=L/2+l(T=T,i j=L/2+l 

Ci Ci = ^L/25'2'/2+l,-f'5'^/2,T ' 

C2 C2 = |tL|5'^|(Si,| ■ ■ ■ Si/2,|)S'^j(SLj ■ ■ ■ 82,72+1, t) , 

and six more terms Cj ® Cj, i — 3, . . . , 8, which are equal to Ci ® Ci and C2 ® C2 with | 
replaced by J, and the hermitian conjugates of these operators. This gives us the inequality 
of the lemma for the lowest eigenvalue of H. The elastic energy terms in the functional £l 
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add up to the same contribution on both sides of the inequahty and therefore can be included 
trivially. Q.E.D. 



Proof of Theorem 1 : 

First we prove that the minimum of £l is attained in a configuration of tj with total fiux 
$ when L = 2 mod 4, and with $ = tt when L = mod 4. Let {tj} be minimizing. Due 
to the invariance of the ground state energy of H under gauge transformations, the energy 
functional Sl defined in (1.3) depends only on {\tj\} and $. Hence, for any set of tj's 
there are t'j such that SL{{tj}) = SL{{t'j}) and t'j^^^ ^'l ^^^^ have the correct 
sign for application of Lemma 7. The two configurations in the right side of (3.1) have fiux 
$ = argument (t^^2^L) which, by assumption, takes the values stated in the theorem. Because 
{tj} is minimizing, (3.1) must be an equality. Therefore the two configurations on the right 
side must be minimizers of £^1,, but then both have the correct fiux, or tt. 

Next we show that in both cases the minimum of Sl is attained in a configuration with 
{\tj\} dimerized. For any configuration {t^^\t^^\t^^^) of real tj's and t*-"^-* of the right signs, 
(3.1) implies that either {t^^\t^'^\t^^^) or (t^^), t('^)) has at least as low an energy and 
has at least as many pairs of identical |tj|'s. In particular this shows that £l is minimized in 
a configuration with \tj\ = \tj+2\ for j = L/2 — 1 and j = L — 1. Because Sl is translation 
invariant the argument above can be repeated to show that there is a minimizing configura- 
tion with \tj\ — \tj+2\ for all j, which is the desired result. Q.E.D. 



Lemma 7 also permits us to solve the fiux phase problem for even rings with U 0. 

8 Corollary. Let Uj = U and let {tj} be a fixed configuration of nonnegative tj 's. Then the 
minimum of Sl, varying over the total flux $ alone, is attained for $ = L = 2 mod 4 
and for ^ — n if L — mod 4. 

Proof of Theorem 4 : 

The Heisenberg antiferromagnet on an arbitrary bipartite lattice is refiection positive [51]. 
This implies, via Lemma 14, the following inequality, analogous to Lemma 7: 

MHueisX{jff,jlr^,4^}) 

where {^fp, Jij^\ "^ij ^} denotes a partition of the coupling constants into three groups: the 
jj^j^ are the couplings on the bonds to the left of any refiection plane of the lattice, the 

(r) 

are the couplings on the bonds intersected by the refiection plane, and the J^j are the cou- 
pling to the right of the plane. Using this inequality Theorem 4 is then proved in the same 
way as Theorem 1, but this time there is no need to distinguish between L = mod 4 and 
L = 2 mod 4 — or even to restrict ourselves to one dimension. Q.E.D. 
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4. ASYMPTOTIC DIMERIZATION FOR L = mod 4; 
PROOFS OF THEOREMS 2 AND 3 

Even though finite rings of length L = mod 4 can do more comphcated things than dimerize 
(see Section 5 for a discussion of counterexamples), they do dimerize asymptotically. To be 
precise, for a fixed /, satisfying conditions i) and ii) of (1.6), for which the infinite ring has a 
unique energy minimizing configuration among the configurations of period 2, any sequence 
of minimizers for finite rings has to approach that dimerized configuration uniformily when 
L oo, i.e., any nearest pair of |tj|'s converges to the pair occurring in the infinite volume 
dimerized configuration. The aim of this section is to prove this statement, which is Theorem 3 
i) . At the same time we will also supply the proofs of Theorem 2 and Theorem 3 ii) . 

In contrast to the previous section, here we will always consider minimization of the 
energy over arbitrary configurations of \tj\, but with constant flux For this purpose it is 
convenient to introduce an explicit notation for the energy functionals at fixed 

Here, and in the rest of this section, it is assumed that all tj are nonnegative. Note that the 
tj 's now play the role of the |tj|'s. 

For any x and y let d{x, y) denote the configuration with tj = x for j odd and tj = y for j 
even, and define the energy density e of a period 2 configuration d{x, y) as in (1.7). Note that 
e is independent of In order to remove the degeneracy due to the obvious symmetry under 
cyclic permutations of the t^'s, we will henceforth adopt the convention that the maximum 
value of all t's is attained by ti. 

A first ingredient in the proof of Theorem 3 is the following energy estimate, which by 
itself implies Theorem 3 ii). Let Ef^ denote the value of in a minimizing configuration 
y Under the general condition (1.6) ii) on / we can find a constant K, depending only on 
/, such that \t^ j\ < K. 

9 Lemma. Assuming only condition ii) on f in (1.6), we have for all — tt < $ < tt 

(\El-El\<2^'^K/L z/L = 2mod4 

\\El- El\<2{'n -^fK/L i/L = 0mod4 

and 

El -Et{tl)\< Att'^K/L ifL = mod 4 
Et - ^i{tt)\ < 8tt^K/L for all ^ and all even L. 

Because Sf^ attains its minimum Ef^ in a dimerized configuration (L = mod 4) , this 
lemma implies, for $ = 0, that even if the true minimum is not dimerized, its energy differs 
from the best dimerized configuration by 0(1/L), which is 2 orders down from El = 0(L). 
Proof : First note that H{{tj}) has real matrix elements and therefore has a ground state 
ip^ which is real, and hence such that {ip^ \ cj^.]^ ^ \ ip^) is also real. Using a gauge 
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transformation that makes H({exp{i7r/ L)tj}) real one can also see that H{{exp{i7r/L)tj}) 
has a ground state ipl such that 

with real Rj,a- 

Let {t^L j} be any minimizing configuration of £^£, and let ip^ be a real ground state of 
H{{t^j}). Then, by the variational principle, 

Et < (Vi I H{{eMi^/L)tlj}) I Vi) + 

j 

= {^L I H{{eMi^/L)tl,}) I Vi) + E/(*?) 

The sum of the first two terms in the right side of the equality above is E^. Using the reality 
of tj)^^ the last term can be estimated by 4|1 — cos($/L)| 'Y^ - 1^5. jl- cosx > 1 — /2, 
this proves the first estimate of the Lemma. The other inequalities are derived in the same 
way. Q.E.D. 



The next lemma is an application of a version of the Abstract Chessboard Estimate [52] 
which is given in the Appendix. 



10 Lemma. For any configuration of tj > 0, with d{tj,tj+i) described in (1.7), 
1 ^ 

JJ2^Lid{tJ,t,+^))<£l{{tJ}) z/L = 2mod4 , 



jJ2^L{d{tj,tj+i))<Sl{{tj}) i/L = 0mod4 . 



L 



Proof : Let $ = if L = 2 mod 4 and $ = tt if L = mod 4. We apply Theorem 15 (with 
D = M^) to the following F: 

F(ai, . . . , ol) = exp(-£:f (ti, . . . , t^)) , 
where aj = (tj.Sj) G IR^. The involution R is defined by R{t,s) = (s,t). The domain V 
is taken to be the set of all (ai, 02, . . . , a^) of the form {{ti, ^2), (^2, ^3), ■ ■ ■ ,{tL, ti))- The 
necessary invariance properties of T> are trivially satisfied. The cylicity of F follows from the 
invariance of under cyclic permutations of the tj . The inequality of Lemma 7 translates 
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into Schwarz's inequality for F on V. This becomes obvious with the definitions: 

a\ = iti,t2) bi = {tL,ti) 

a2 = {t2-,h) &2 = (*l-i,^l) 



t*L/2 — ^L/2+l) ^>L/2 — (^L/2+1: ^L/2+2) • 

Indeed, one then has: 

exp(-^f(ti, . . . , tL)) = F{ai, a2, . . . , 0^/2, ^>l/2, • • • , 

exp{-Sf{ti, . . . , tL/2,tL/2+l, tL/2, • • • , ^2)) = F{ai,a2, . . . , aL/2, RO'L/2, • • • , -Rtti) 
exp(-£^|'(tl, tL, • • • , tL/2+2, tL/2+1, • • • , tL)) = F{Rbi, Rb2, . . . , RbL/2: bL/2: ■ • • , 
exp(-£:f (tL, ^L-l, ■ • • , tL/2+1, ^L/2+2, • • • , ^l)) = ^(^1, ^2, • • • , &L/2, ^^L/2, • • • , ^^l) 

The imphcation of Theorem 15, in terms of then reads: 

^L (^) ^ ^ ^* ' ' *J ' ' • • •) ' 

i=i 

which is the desired inequahty. Q.E.D. 



11 Lemma. For the given f let K be a constant such that Itfjl — ^> f^''^ ^- ^^''^ 
L = mod 4, let {t^j} be any minimizing configuration of Then 

md{tl„tl^^,))-El\<4K^ , 

for allj = 1,2, ...,L. 

Proof : As L = mod 4 we can apply Lemma 10 to and with tj =t^j- We then have 

^i2^Lid{tl,,tLj+,))<Sl{tl) . 
By Lemma 9 the right side is bounded by El + Att'^K/L, and hence 



\Y.{^im,,tl,^,))-El]< 



L 

As each of the terms in the sum is non-negative, the statement of the lemma follows. Q.E.D. 



12 Lemma, i) The following limit defines a continuous function e: 

lim Sf{d{x,y)) = e(x,y) 

L — >oo 

The convergence is uniform on compact subsets of IR^ . 

ii) If e attains its minimum at {xQ^yo) and (yo,^o) o-'^^d nowhere else, then there exists a 
sequences xl — xq and yL — J/o such that for all L = mod 4, d{xL,yL,) is a minimizing 
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configuration for S^, and such that 

lim Y£l{d{xL,yL)) = e{xo,yo) 

i— »oo L/ 



Proof : The arguments involved in proving this lemma are fairly standard. We therefore 
only give a sketch of the proof 

i) We first show that there exists a concave function e{x^ y), satisfying e{x, y) < 2|a;|+2|j/| + |t/|, 
and such that 

e{x,y) = Yim. \\Q{HL{d{x,y))) . 



For brevity we put Gl — ^o{HL{d{x,y))). By the variational principle and the fact that 
||cj>|| = 1 one immediately gets that Gl is "almost subadditive": 

Gl,+l^ < Gl, + Gl, + 6(|a;| + \y\) 
and therefore GL + 6(|a;| + |j/|)is subadditive. This sequence is also almost monotone: 

Li > L2 ^ Gli < Gl2 + Constant 
because Gl < 0. It follows that the limit of Gl/L exists and as a limit of concave functions 
the limit is also concave. The bounds are trivial but they imply continuity and therefore the 
convergence is uniform on compacts. From Lemma 9 it is obvious that e does not depend on 

ii) From the conditions on / (1.6) one has an apriori bound K on the minimizing tj: I^JlI — ^■ 
By Theorem 1, always has a dimerized minimizer d{xL,yL), and by compactness there 
always exists a subsequence (xL^^yLk) converging to By the continuity of e'^ and 

the uniform convergence on compacts 

■^£LM(^Lk,yLk))-eix^,y^) ^0 . 

As the d{xLf,,yLk) are minimizers it follows that e(a;*,|/*) = e(a;o,yo) and hence, by assump- 
tion, (a;*,!/*) = (a;o,2/o) or (a;*,?/*) = {yoiXo). Without loss of generality we can assume 
> yQ. Consider then the new sequence d{xL,yL) of minimizers with {xL,yL) = {xLiyL) 
or (yLjXL), and such that XL>yL for all L. Then all convergent subsequences (xL^^yLk) of 
{xL,yL) converge to the same limit {xo,yo) and hence the sequence is convergent. Q.E.D. 



Proof of Theorem 3 i): 

Let d{xL,yL) be a sequence of minimizing configurations for SJ^ and such that xl and yL 
converge to xq and yo respectively as L — > 00. The existence of such a sequence is guaranteed 
by Lemma 12. 

First we estimate e(t° ^lj+i) ^ minimizing configuration t° of and any 
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< 



+ ^ \^Iiditlj,tlj+^)) - £lidixL,yL))\ + ^ \ElidixL,yL)) - e(a;o,yo)| • 

All three terms on the right side vanish as L — > oo, uniformly in j: the first term by the 
uniform convergence on compact sets given by Lemma 12 i), the second by Lemma 11, and 
the third term by Lemma 12 ii). 

The theorem now follows from the uniqueness of the dimerized minimum and compact- 
ness. Indeed, suppose one can find arbitrarily large Lj. such that 

\tl,,j-x\ + \tl,,j+i-y\>^ and \tl^,j-y\ + \tl^,j+i-x\>s , 

for some e > 0. By compactness and the continuity of e there would then be a limit point 
{x', y') of {t'j^^ J, t^^ j+i)' ^'-'^ which e{x', y') = e(xo, yo). As (x', y') as at least distance e away 
from [xq, 2/0 ) this would imply the existence of another minimizer for e which is ruled out by 
assumption. Q.E.D. 

Proof of Theorem 2 : 

Let $ = if L = 2 mod 4 and $ = tt if L = mod 4. Denote by d{x, y) the unique dimerized 
minimizer of S^. Let {tf} be any minimzer of S^. As attains its minimum in {tf }, the 
inequality of Lemma 10 is an equality for tj = tj . Therefore, all configurations on the left 
side of (4.1) must be minimizers. These configurations are all dimerized and therefore (by 
our uniqueness assumption) must all be equal to d{x,y). Q.E.D. 
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5. COUNTEREXAMPLES: RINGS OF LENGTH L = mod 4 

HAVING MINIMIZERS AT ZERO FLUX THAT ARE NOT PERIOD TWO 

We will only discuss counterexamples with U = 0. It is obvious, then, that counterexamples 
also exist with small non-vanishing U. However, the presence of a not too large on-site 
repulsion tends to enhance the dimerization, while nearest neighbour repulsions can reduce 
it. We have not investigated the effects of interactions on the new instabilities discussed in this 
section. Our result for the spin-Peierls problem. Theorem 4, which can be interpreted as the 
U oo limit [26], seems to indicate that, when U is very large, the translation symmetry 
can break down from period one to at most period two. This seems natural because the 
electrons are more mobile when U is small and, under this condition, other deformations 
of the lattice can become favorable if the flux is not optimal. We interpret distortions that 
break the periodicity two as an attempt of the system to mimic the energy lowering effect of 
a nonvanishing flux. This is also a possible interpretation of the lattice distortions found in 
two dimensions [45,46] . 

Kennedy and Lieb [18] gave an example of a function / (involving a quartic term) which 
gives rise to a minimum for the model on a square (L — 4) which does not have periodicity two. 
Longuet-Higgins and Salem pointed out the possibility of other than period two instabilities 
a long time ago [15], but they expected these to be unimportant and thought the period two 
instability would always dominate. They only considered the case L = 2 mod 4 and therefore 
did not notice that other instabilities can, in fact, occur when L — mod 4. 

Our aim in this section is twofold. First we present a family of functions / that lead 
to non-period- two minimizers of the zero-flux problem for the ring of eight sites. Then we 
discuss what instabilities other than period two occur. From the discussion it will become 
clear how to construct many more examples of non-period-two minima. The properties of 
the new instabilities, which we study partly only numerically, support our conjectures stated 
in the introduction. 

Counterexample for L = 8. 

For simplicity we will present just one family of functions / chosen to yield simple values 
for the minimizing dimerized configuration. It is easy to generalize this example in several 
directions. 

Define a function f{t), for t > 0, as follows (see Figure 2 for the graph of /.) 

r ± + l{t-l) + ^{l-tr ff0<4-a 

m - (5.1) 

(± + l(t-l) + ^il-t)'+^it-4 + a)^ ift>4-a . 

Here a can be any number strictly between and 4 — 5/a/2 = .4645. 

Consider first configurations tj of the form t2j+i = t2j = y. We always assume 
X > y > 0. The ground state energy of eight electrons on a ring of eight sites described by 
the Hamiltonian (1.1) with U = is then 4e(x, j/) where e{x,y) is given by 

e{x, y) = —X — \/ -\- %p- . 
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For the 8-ring with / given in (5.1), the optimum period two configuration is then determined 
by the values of xq > yo that minimize 



h{x, y) = e{x, y) + f{x) + f{y) 



It is trivial to see that xq > 0, and as h is difi^erentiable away from x = and h{x, y) — > +oo 
as a; or J/ —> +oo, {xo,yo) must therefore be a critical point of h{x,y) solving the equations 



A straightforward calculation shows that {x, y) = (4, 3) is the unique minimum. 

Next we study the stability of the relative minimum {xq, yo) under a class of perturbations 
which break the period two symmetry. For future use we do this for arbitrary L = mod 4, 
L > 8. More specifically we consider configurations of the form t2j+i = t2j = y + ^Vj+i, 
j = 1, . . . ,N mod N, with L = 2N. As we do not want to introduce a flux, the Vj must be 
real. One can use analytic perturbation theory to compute the energy as a function of x, y 
and s. Up to corrections of 0{e^) one finds 



X 



1+ 



y 




f\y) 



Kin. Energy 



2 ^ An - £ 




n=l 



N / 



|V"(0)P(y +XCOs(Q!n))2 Zlm7^2n mod at 





(5.2) 



where 



An = -x/a;^ + j/^ + 2a;ycos(Q!n), a = 2tt/N 



and 



n,m — 



y^ + x"^ cos{a{n + m) /2))^ + 2xy cos{a{n + m) /2) cos(Q;(n — m) /2) 



V is the Fourier transform of Vf V{q) — {l/N) exp{27rijq/N)Vj. Note that the coeffient 
of £^ is diagonal in V. Up to terms 0(£^), the elastic energy is 
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In a critical point (x, y) in the set of period two configurations, the terms proportional to 
e in the electron energy will cancel the ones in the elastic energy. A perturbation V will 
then lower the energy for small e if the total (electron plus elastic energy) coefficient of 
is negative, which means that the optimal dimerized configuration is not the true minimum. 
Call the coefficient of the contribution to the electron energy —A{x,y). Note that this 
coefficient is always negative. There is an instability whenever 

N 

A{xo,yo)>Yf"{yo) . (5.3) 

We now return to the model on the 8-ring and show that indeed the energy can be 
lowered by perturbing the dimerized minimum (xq, yo) found above. For the 8-ring (N = 4) 
there are four independent perturbations Vi to consider: 

F = (f,f,f,f), F = A/2(f,0,~f,0) , 
V = (1,-1, 1,-1), V = ^(0,1, 0,-1) . 

The first, V = (1, 1, 1, 1), respects the period 2 and in the optimum dimerized configuration 
it cannot lower the energy. The second and the fourth have the same energy because they 
are translates of each other. We have normalized the V^s such that Yli "^"i — 4 (in general 
we will use the normalization 'Ylii^i — ^ ■) Let — Aq, — Ai, — A2 be the coefficients of 
for the first, second and third choice of V respectively. It is convenient to factor out a; as a 
scale factor and to use the ratio r = yjx as the relevant variable, < r < 1. The explicit 
expressions for the A's are: 

- 2 

^^0- (l+r2)3/2 

2 

xAi = 4 



J- -r 

xAo = 2 



In our example both Ai and A2 satisfy the criterion for instability (5.3): 

|Ai(4,3)| = ^>2r(3) = | |A2(4,3)| = ^>2r(3) = ^ . 

Of course, Ao(4,3) does not satisfy (5.3), indicating that (4,3) is a stable relative minimum 
within the set of period two configurations (indeed Ao(4, 3) = 32/125 < 2/5). We conclude 
that for the choice of / defined in (5.1) the energy minimizing configuration of the 8-ring is 
not period two. 

Other than period-two instabilities for arbitrary L = mod 4. 

The behavior of the coefficients of in (5.2) for large N gives support to our conjectures 1 
and 2 stated in the introduction. 
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Let -Afc denote the coefficient ofs^ for V = V^'^\q), where V^''\q) = {Sq^k+Sq,N-k)/ V^, 
k ^ 0,N/2, V(o)(g) = Sq,o, V^^/^\q) = Sq^N/2- The perturbations are normalized such that 
their contributions to the elastic energy is identical. Aq and Ai are given by: 

X smlan) 



^ An Xn+1 J xy{cos{an) - (x>s{a{n+l))) 



The stability of a relative minimum {xotUq) within the set of period two configurations re- 
quires that Aq < {N/2)f"{yQ). When Aq < (A^/2)/"(|/o) < Ai such a relative minimum is 
unstable against the perturbation eV^^\ which has wavelength = L. Numerical investigation 
of these expressions reveals the following: for r = y/x < 1 but close enough to 1 (depending 
on N), one has Ai > Aq. This is illustrated in Figure 3. Note that the perturbations with 
k > 1 are all stable. We showed in the previous subsection that for L = 8 this situation can 
be realized with a function / satisfying the conditions i and ii) stated in the introduction. It 
is reasonable to expect that this can be done for all L that are a multiple of 4, but the larger 
L becomes the more stringent are the restrictions on /. In particular, if one fixes / we expect 
that for large enough systems the long wavelength instability will not occur and the system 
will remain in a period two configuration. 



A. Appendix 

1) Discontinuous and non-convex / 

We have assumed throughout this paper that / was a continuous function satisfying some 
simple conditions as stated in the introduction (1.6). Whether or not / is continuous we 
can associate a function / with / having the following properties: (a) / is convex (i.e., 
/(At + (1 - X)t') < Xf{t) + (1 - X)J{t') for all < A < 1 and aU t, f > 0; (b) J{t) = sup{g{t) : 
g is convex and g{t') < fit') for all t' > 0}. This function / is called the convex hull of / 
and is automatically continuous for t > 0. The following lemma shows that all we really need 
is that / (not /) satisfy conditions i) and ii) of (1.6). 

13 Lemma. Let {tf} be a minimizing configuration of£^. Then for all j, f{tj) = f{tf), 
where f is the convex hull of f. 

Proof : Let j be arbitrary and fixed. As a function of tj = t the energy can be written as 

£{t) = E{t) + f{t) 

= {E{t) + at) + {fit) - at) = Kit) + git) . 
We choose a = (/(t2) — /(^i))/(^2 — ^i) where 

ti = sup{s<tf|/(s)^7(s)}, t2 = inf{s>t*|/(s)^7(s)} . 



24 



As t appears linearly in the Hamiltonian H, K{t) = E{t) + at is a concave function of t (i.e., 
—E{t) — at is convex). By construction we also have that g{t^) > g{ti) = 5'(^2)- Hence, with 
< q; < 1 such that tj = ati + (1 — a)t2 one has 

£{tf) > aS{ti) + (1 - a)8{t2) 
and the inequality is strict if git'j) = f{tf) — f{tf) ^ 0. Therefore the minimium can be 
attained only at tj where /(tj) = jitj)- Q.E.D. 



2) The DLS Lemma for ground states 

The following lemma is a ground state version of Lemma 4.1 in [51]. An example of this 
ground state version is given in [53]. Note that the matrices A and B need not have real 
matrix elements. It is obvious how to generalize to reflections that involve an additional 
unitary transformation. When needed one can in fact apply the lemma as stated below after 
performing a unitary transformation on one half of the system. In this paper we applied the 
lemma in two different situations. In the proof of Lemma 5 the two copies of the Hilbert 
space 7i correspond to spin | and spin J,. In Lemma 7 the two spaces refer to the right and 
left halves of the ring. 

14 Lemma. Let A, Ci, . . . , Cn be a collection of d x d complex matrices (n could he 
infinite) with the following properties: A and B are Hermitian (i.e., A = A^ , B = B"^), for 
all i, Ci is real and Ci ® Ci is symmetric (as a d"^ x d'^ matrix). Let Xo{A, B) denote the 
lowest eigenvalue of the matrix 

T{A,B) = A®l-^Ci®Ci + \®B . (A.l) 

i 

Then 

Ao(A B) > i (Ao(A, A*) + Xo{B\ B)) , (A2) 

where A* denotes the matrix obtained from A by complex conjugation of the matrix elements. 
In particular 

Xo{A, B) > min (Ao(A, A*), Xo{B*, B)) . 

The inequality (A. 2) is strict if for some i 

(O I ® I I O) ^ (O I I ® I O) , {A.3) 
where Q is any eigenvector for the eigenvalue Xo{A,B) (this is not the only case where the 
inequality is strict.) In particular all ground states Q ofT{A, A*) satisfy 

{n\c^(g)l\n) = {n\l(g)C^\n) . 

Proof : Let {ej} be a orthonormal basis of in which the Ci have real matrix elements. 
Define V to be the antilinear map leaving the Cj invariant: Vcj = Cj. Then VAV = A* etc. 
Let O, be an eigenvector of the smallest eigenvalue of the operator (A.l) and denote by 
its coefficients in the basis (8) ej : 

Q = MijCi ® Cj 
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We can consider (My) as a matrix and use the eigenvectors {ua} of M^M and {va} of MM^ 
to define two new bases: 

<Pa = ^{^j I Ua)ej , Ipa = ^{^j I ^a)*ej . 

j 3 

It is then straigthforward to check that 

a 

We can obviously assume that Aq, > by absorbing a phase in the definition of ipa, and by 
omitting any vanishing terms in the sum. 

Using the fact that {^Pa} and {V'/s} are orthonormal bases the energy of Q, can be expressed 
as 

(O I T{A,B) I O) = ^ a2 ((v?, I a I + (V'a I S I V'a)) 

a 

i af3 

This energy can directly be compared with the energy of T{A, A*) and T{B*,B) in the states 
Qi and Q2 defined by 

a 
a 

which have the same norm: = = ||^2p — Yla "^a- ^^^^ average of the energies of 

Oi and O2 is 

I T{A,A*) I Oi) + (O2 I T{B\B) \ O2)) 

= 1 5^ (((^^ I A I (^«) + {Vifa I A*F I (^«) + {VlPa I S*F I V'a) + (V'a I ^ | V'a)) 
a 

- A«A;3 ((V'a I C', I ^p){V^a I C'.y I (/p^) + (V'a I Q I i^p){Vi;^ I I ^A^)) . 

i aP 

{AA) 

Using (y f! \ ip) — {lp \ V"^ \ ipY ^ for any two vectors (/?, '0, the hermiticity of A and -B and the 
reality of the it is straightforward to verify: 

{V(pa I A*V I (/?a) = {iPcc\A \ ipa) 
{Vipc I CiF I ipp) = {ipa I Ci I V?/?)* , 

etc. As Aq, > we can then apply 2\uv\ < \u\^ + \v\^ to each term in the double sum of (A. 4) 
to obtain 

{VL I T{A, B)\VL)>^ {{VL I T{A, A*) \ VL) + {VL\ T{B*, B) \ n)) 
from which (A. 2) directly follows by the variational principle. 

The inequality is strict if for some i and some a,P, {ipa \ Ci \ ipp) 7^ (-^q, | Ci \ if^p), 
for which (A. 3) is a sufficient condition. When A = B*, (A. 2) is an equality and so 
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{(Pa. I Cj I (fp) ^ {tpa I Ci I tpp) for all i, a, (3 in that case. Q.E.D. 

From the proof of the lemma it also follows that any Hamiltonian of the form T[A, A*)^ 
has a ground state VL satisfying 

(fi I X X* I f2) > . 

for all observables X orvl-i. 

The same state therefore also satisfies 

{Q\X = {n\l®X* \ Q) . 

In particular if [X ® I, T(A, A*)] = 0, and X = X"^ , there will be a ground state O such that 
X ® Ifl — 1 ® X*fl — xfl for some eigenvalue x of X . (This also follows from application of 
(A. 2) with A replaced by ^ + c{X — xi)"^ and B hy A* + c{X* — X2)^, where xi and X2 are 
the eigenvalues of X (8) I and I (8) X* in some ground state of T{A, A*). 

3) The Abstract Chessboard Estimate([52]) 

15 Theorem. Let D be a set, let R : D ^ D be an involution on D (R"^ = identity), and 
let F : V (D be a complex valued function on a domain V C D^'^^ that has the following 
invariance properties: 

(ai, . . . , a2iv) eV^ (a2, • • • , a2N, ai) E V , 
(ai, . . . ,aN,bN, ■ ■ ■ ,bi) eV ^ (ai, ...,aN, RaN, ■ ■ ■,Rai), (6i, . . . , ^at, RbN, ■ ..Rbi) e V 
Assume F obeys 

F{ai, . . . , a2jv) = F{a2, a2N, oi) 

and 

\Fai, . . . ,aN,bN, ■ ■ ■ , bi)f < F{ai, ...,aN, RaN, Rai)F{bi, ...,bN, RbN, Rbi) . 
Then \\a\\ = \F{a, Ra, A, . . . , Ra)\^/^^ satisfies 

2N 

\F{ai, a2N)\ < Y\_ W'^iW ■ 

The proof of the Theorem is exactly the same as in [52] Theorem 4.1. In that reference 
it is also assumed that D is a real vector space and that F is defined on all of D^^ . Under 
these circumstances the inequality of the theorem implies that || • || is a seminorm on D, but 
we do not use this property in this paper. 
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FIGURE CAPTIONS 



Fig. 1. A ring containing an even number of sites, labeled 1, . . . , L. The hopping matrix 
element tj is associated with the bond {j,j + 1}. Reflection through the plane indicated by 
the dashed line interchanges the role of the hoppings labeled t^^^ and the ones labeled t^'^^ 
(see Lemma 7). 

Fig. 2. Graph of the function f{t) deflned in (5.1) for a = .46. Recall that large t corresponds 
to small separation w and vice versa. 

Fig. 3. Graph of the coefficients xAj~ as a function of the ratio r — y/x for L = 20 and 
k = 0, ...,5. An instability can occur when Afc(r) > Ao(r). Note that this condition is 
satisfied only in a small region of r not much smaller than 1, and only for k = 1, which 
corresponds to an instability of wavelength L in real space. 
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